Investigating quantum transport with an initial value 
representation of the semiclassical propagator 

Christ oph- Marian Goletz, Frank Grossmann 
Institut fiir Theoretische Physik Technische Universitdt 
01062 Dresden, Germany 

Steven Tomsovic 
Department of Physics and Astronomy 

Washington State University 
Pullman, Washington 99164-2814, USA 

Abstract 

Quantized systems whose underlying classical dynamics possess an elaborate mixture of regular 
and chaotic motion can exhibit rather subtle long-time quantum transport phenomena. In a short 
wavelength regime where semiclassical theories are most relevant, such transport phenomena, being 
quintessentially interference based, are difficult to understand with the system's specific long-time 
classical dynamics. Fortunately, semiclassical methods applied to wave packet propagation can 
provide a natural approach to understanding the connections, even though they are known to 
break down progressively as time increases. This is due to the fact that some long-time transport 
properties can be deduced from intermediate-time behavior. Thus, these methods need only retain 
validity and be carried out on much shorter time scales than the transport phenomena themselves 
in order to be valuable. The initial value representation of the semiclassical propagator of Herman 
and Kluk [M. F. Herman and E. Kluk, Chem. Phys. 91, 27 (1984)] is heavily used in a number 
of molecular and atomic physics contexts, and is of interest here. It is known to be increasingly 
challenging to implement as the underlying classical chaos strengthens, and we ask whether it is 
possible to implement it well enough to extract the kind of intermediate time information that 
refiects wave packet localization at long times. Using a system of two coupled quartic oscillators, 
we focus on the localizing effects of transport barriers formed by stable and unstable manifolds in 
the chaotic sea and show that these effects can be captured with the Herman-Kluk propagator. 



PACS numbers: 03.65.Sq, 05.45.-a 



I. INTRODUCTION 



Model Hamiltonians for many molecular and atomic systems can be represented in the 
form of coupled anharmonic oscillators . The dynamics of such systems, generally 

speaking, are not fully integrable, and may be said to be in the near-integrable regime if, 
for the most part, the resonances introduced into the dynamics by anharmonicities are non- 
overlapping Si- On the other hand, if a great enough number of the resonances overlap 
strongly, a well developed domain of chaotic dynamics is generated, which coexists with 
other domains of near-integrable behavior. In this case, the dynamics is in the mixed regime 
and the focus of the investigation presented here. Due to the underlying dynamical diversity, 
quantized versions of such systems may display many remarkable properties and effects: for 
example, dynamical tunneling , chaos-assisted tunneling 0,0, S] and ionization flo| . 
and various forms of transport suppression and eigenstate localization 
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Semiclassical theories, which date back to the early days of quantum mechanics 20|, 



22, 1231, 12^, provide a powerful approach to understanding the quantum dynamics of 
mixed systems, yet are incompletely developed. Although trace formulae, which connect 
quantum densities of states to sums over periodic orbits, exist for integrable 
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281], and chaotic systems 29|, analogous formulae are not known for mixed 



systems. The study of dynamica' 



fully understood; see however [30 
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chaos still remains to be 



39| for progress. Similarly, 



the time-dependent or dynamical basis of various mechanisms of eigenstate localization are 
poorly understood, and there has not been a consensus that such localizing effects are 
captured in a full semiclassical approach 
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42] . The discovery of a sum rule for subtle 



classical orbit action correlations seems to support the studies which find the semiclassical 



theories capable of describing aspects of localization 
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However, it is not necessary to resolve this question to deduce useful information. For 

presence of classical transport barriers 
(for example, cantori) intermediate time scale dynamics are sufficient to deduce aspects of 
long-time dynamical properties such as suppression of transport. It is known that even for 
strongly chaotic systems, non-uniformized implementations of semiclassical dynamics are 
accurate to time scales that lengthen much more in the short wavelength limit than the 
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time for chaos to develop on an ever finer scale [45|, |46j. This is usually expressed in terms 
of a scaling with Planck's constant h. For smooth Hamiltonians, such as strongly coupled- 
anharmonic oscillators, the time scale of accurate semiclassical propagation would lengthen 
algebraically as h^^^^ 47j or better. For integrable systems, the expectation is at least as 
long as h~^. Mixed systems can be made to work at least as well as chaotic systems. The 
time scale for chaos to develop to the course-grained volume represented by a quantum state, 
roughly h'^ (where d is the number of degrees of freedom), only grows logarithmically as r oc 
In h~^, where fi is the Kolmogorov-Sinai entropy (essentially related to positive Lyapunov 
exponents). So although transport suppression (or localization) time scales would roughly 
be much, much longer, oc h~'^, the dynamics between the logarithmic and relevant algebraic 
time scales does provide a window wherein interference is dominating the dynamics and with 
which there is at least an opportunity to deduce something about long-time properties such 
as quantum transport suppression or eigenstate localization. 

In the phase space region dominated by chaotic dynamics, it is quite common to find 
partial barriers that limit transport connections between subregions of the chaotic domain. 
The barriers can be constructed as surfaces having the property that locally a minimum flow 
of phase points cross from one side to the other per unit time; i.e. a minimized flux. These 



bounding surfaces can take the form of so-called cantori 48|, |49| or intersections of stable 
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52|, both of which in two-degrees-of- 



and unstable manifolds of short periodic orbits 
freedom are associated with turnstiles that determine quantitatively how impeded the flow 
is. They introduce additional time scales into the dynamics. The semiclassical theory must 
be developed and valid long enough to be able to deduce these time scales. 

Semiclassical theories can be based on stationary phase approximations applied to various 
representations of the Feynman path integral. In atomic, chemical, and molecular physics 
applications, the so-called initial value methods 
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551 ] are heavil y re lied upon, more 



so than the boundary value methods, such as in the work of van Vleck 2J] . It is of interest 
to understand better the limits and strengths of these initial value methods. Here, we focus 
on the implementation derived by Herman and Kluk (HK) 5J] which was shown to be 



very efficient. T 
approach 
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le dynamics of molecules and atoms can be treated successfully with this 
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591], although the more unstable the dynamics in any chaotic domain. 



the more challenging the work. Here, we apply it to the semiclassical propagation of wave 
packets in chaotic regions of a two-dimensional mixed phase space quantum system. It 
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turns out that the HK method is suitable for propagation times where locahzation due to 
transport barriers becomes manifest. In contrast to previous works of Yoshida et al. [4^ 
and Maitra 42|], we focus on obtaining the critical intermediate-time dynamical information 
of a mixed system in a parameter regime where the chaotic domain is strongly unstable, yet 
nevertheless retains significant partial transport barriers. 

This paper is structured as follows: in Sect. [TTlwe review using the propagation of localized 
wave packets as a means of investigating the suppression of quantum transport due to 
partial classical transport barriers in phase space. Intermediate propagation time dynamics 
is connected to suppression of quantum transport seen at long times with a combination of 
semiclassical, random matrix ensemble, and master equation approaches. This is followed by 
a description of the semiclassical initial value representation and purely classical methods in 
Sect, mil The coupled quartic oscillator model system is introduced in Sect. [IVl In Sect. IVl 
our numerical results are presented and followed with a brief discussion and conclusion. 

II. INHIBITING QUANTUM TRANSPORT 

Because it has not received a great deal of attention and it is central to our work, we review 
in this section previous work demonstrating the connection between short- and intermediate- 
time propagating wave packet behaviors and long-time, inhibited quantum transport due 



to partial classical transport barriers 18|, iGOj. Of course, the long-time suppression of the 



spreading of a wave packet implies a localization of the eigenfunctions of the system and vice 
versa so that these two topics are manifestations of the same physics and discussion of one has 
its translation for the other. The weakest forms of localization related to the chaotic zones 



in phase space, i.e. coherent backscattering/weak localization 6l| and eigenstate scarring 



by short, unstable periodic orbits [62], are specifically excluded from the approach presented 
in this section through the assumptions regarding local mixing and the absence of certain 
correlations. However, they have been discussed at length elsewhere, would increase the 
localization discussed below, but are not of primary importance here. 



4 



A. Transport barriers' connection to quantum states 



Consider a system with a mixed classical dynamics such that on any energy surface of 
interest, there exist transport barriers which can be used to separate the full chaotic phase 
space part into several subregions. Following Bohigas et al. 18|, a first critical concept in 
relating these classical transport barriers to quantum behaviors is defining state vector sub- 
spaces associated with each separate region of the available chaotic phase space in a narrow 
energy range of interest. It is with respect to any basis, which respects these subregions, that 
gives one a foundation for discussing eigenfunction localization in the dynamical context. 
There are a few parameters containing the most critical information, i.e. the j*^ subregion's 
phase space volume Vj in a narrow (Thouless) energy window AE 63|, vector subspace di- 
mensionality Nj ^ Vj/h*^, relative fraction of phase space volume fj = Vj/V, and a shortest 
mixing time scale oc In (Vj/h'^). Assuming all the subregions' volumes are the same 
order of magnitude, is the time for the chaotic mixing in the dynamics to take a cell of 
volume h'^ to mix into all the cells within one of these chaotic subregions. A transport barrier 
is relevant if its flux (defined on a fixed energy surface) crossing from the j*'* region into the 
f^th j-egion, is small enough that Vj/cpjk » TmAE; generally (j)jk is energy dependent, 
i.e. depends on the local energy surface in question. From this, one sees already that in 
the short wavelength limit (summarized as h ^ 0), any single barrier loses its ability to be 
responsible for localizing effects with respect to the defined subregions because in this limit 
eventually TmAE » Vj/(j)jk- However, this limit is approached extremely slowly due to the 
logarithmic dependence. In a broad range of applications, transport barrier time scales are 
found to be quite relevant. In the case of diffusive dynamics, an infinity of barriers may be 



present, each o 
localization 64 



which is nearly negligible. This is intuitively a dynamical basis for Anderson 



65|. 



It is the mixing time scale that underlies the possibility of a combined semiclassical 
and statistical approach. In the vector subspace associated to a single chaotic subregion, 
quantum states that are approximately stationary (up to times of Vj/[(j)jkAE]) have the 
statistical properties implied by random matrix theory according to the Bohigas- Giannoni- 
Schmit conjecture 661] (random matrix theory provides a proper way to define ensembles 
with which to define statistical measures). Crudely, these states can be thought of as "un- 
perturbed" eigenstates of a system in which the transport fluxes crossing from one region 
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to another are vanishing. Thus, the quantum evolution operator is described by a matrix 
that divides into n x n blocks, where n is the number of subregions. Each diagonal block 
would have the properties of the unitary random matrix ensemble if all the fluxes (pjk van- 
ished. In that case, except for tunneling tails, any initial wave packet essentially within 
one subregion would remain within its respective subspace for all propagation times, and 
the system's Hamiltonian expressed in this representation would also have diagonal blocks; 
its spectral statistical properties would mimic a weighted random superposition of random 
matrix spectra. For no n- vanishing (pj^, the evolution operator and the Hamiltonian would 



eak 

■ 3, 



have matrix elements connecting the diagonal blocks, and thus a wave packet would 
out at least partially from the subspace to which it initially belonged. In French et al. 
it was stressed that there existed a universal, dimensionless transition parameter A for sym- 
metry breaking problems (here a dynamical symmetry is being broken) defined as the mean 
square admixing matrix element of the Hamiltonian divided by the mean level spacing 
squared D^. They emphasized that although it is straightforward to identify the transition 
parameter through perturbation theory, it had a far more general meaning and contained 
the necessary information well beyond the validity of perturbation theory. For this prob- 
lem, Bohigas et al. [l^ showed semiclassically that the quantum transition parameter could 
be expressed in terms of the classical fluxes and relative fractions of phase space volumes 
through the relation 



A _ '"jk _ (Pjk X 



i 

and it determines the degree to which dynamical localization is present in the system. For 
very small values of the transport flux relative to /i, presumably this formula should be 
modified to incorporate a tunneling contribution. However, to our knowledge the tunneling 
correction is still an unsolved problem. 

It is worth noting here a relationship between a quantum transition probability and the 
transition parameter. In a perturbative regime, for short times (but still longer than the 
logarithmic mixing time scales), and a single barrier, the relationship is given by [3] 

- - - 27r 

{u\t)p^u{t)Pu)E = j-^^hh^^kt. (2) 

where p{E) is the local mean level density, lJ{t) is the unitary evolution propagator, and 
the {Pj, Pfc} are projectors into the subspaces defined above; as indicated, the expectation 
value is taken locally in energy. The quantity on the l.h.s. above measures how much of the 
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k*^ subspace leaks into the j*^ subspace under evolution. This can be envisioned as a kind 
of Fermi-Golden-Rule relation for chaotic systems. 

B. A measure of quantum transport suppression 

In order to investigate quantum transport using the dynamics of wavef unctions, it is 
convenient to consider wave packets \E'j which reside within some particular vector subspace 
associated with a (say the j*^) phase space region. The propagations of these initially 
"localized" wave packets are then studied by observing their projections onto the 

various other localized wave packets as a function of time. A particularly convenient measure 
is given by 

T+AT 

^^■'^ = aP^^oo ^ / K^^l^Mfdt, (3) 

T 

which is the long time average of the absolute square of "^jit) projected onto some "^k- 
The notation (■) indicates the appropriate norm (overlap). Assuming the partial barriers 
are effective, ^E'j(t) will spread with different mean intensities into each region, which will 
saturate at some point and not increase any further no matter how long a propagation time 
considered. Statistically speaking, it is possible to deduce the extent of such a localization 
property from the Vjk with one extra ingredient, i.e. the smoothed local spectra of (and 

S,{E) = J2\m3)fSe{E-E^) (4) 
1=1 

where the are the previously mentioned approximately stationary states of the j*'' 

region (unperturbed eigenstates) and e indicates that the (5-functions are smeared sufficiently 
to smooth out oscillations locally in energy. For a Gaussian wave packet not-too-close to 
a turning point, there is an approximate analytic expression for the envelope of the local 



spectrum given by Fourier transforming the initial decay of the autocorrelation function 53| , 
which gives. 



SliE) ^ ^,„„ „ exp ( 'J ) . (5) 



Here 7 is a width parameter, po the mean momentum, and Eq the mean energy. If the wave 
packet is close to a turning point, a more accurate expression should be used. 
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It was shown 



that 



jk 



A,fc(A) 

fkPc 

A,fc(A) 



dE Sj{E)Sk{E) 



(6) 



where pc is the locally smoothed density of states for the full chaotic region. All the infor- 
mation about the extent of localization is contained in the function Ajfc(A) whose meaning 
is given by the projection into the k^^ subspace of a propagated approximately stationary 
state from the j*'' region 

T+AT 



lim — — 

AT^oo AT 



dt 



Pk\m) 



(7) 



In a first order perturbative regime (i.e. small transport fluxes regime), it can be derived 
analytically with the result 18|] 



A,fc(A) = hV^^k (j ^ k) 

A,-,(A) = l-J2fkV^k 

k 

and in the opposite extreme of very strong mixing, it's value saturates at 

lim AjkiA) fk . 

A— >oo 



(8) 



(9) 



The closer Ajk{A) is to the value fk, the less localization in the eigenstates or the less 
suppression of quantum transport. 

The perturbation result for Ajfc(A) can be generalized using a Master equation approach 
developed in Smilansky et al 60|]. It relies on the construction of a Markov matrix whose 
eigenvalues and eigenvectors contain the all the information about transport connections and 
volumes of the subregions. They define a classical probability vector P with n components 
that at any moment in time gives the probability of being found in one of the n subregions. 
Defining the transition probability matrix A, whose elements are given by 



A 



jk 



<Pjk 



(10) 



fc=i -^J 



gives the probability vector at any time 

P{t) = exp(-At)P(0) (11) 

Since total probability is preserved, there is always one vanishing eigenvalue of A whose 
corresponding eigenvector represents the state of complete mixing or uniform probability of 
being found anywhere in the volume. The rest of the eigenvalues and eigenvectors describe 
the decay of any initial probability towards the state of complete mixing. The smallest 
eigenvalue bounds the time it takes for this to occur. The quantum spectral fluctuations 
were related to the trace of the operator exp(— At) using a semiclassical theory. A similar 
calculation could be done for the Ajk{A). 



III. THE SEMICLASSICAL PROPAGATOR IN INITIAL VALUE REPRESENTA- 
TION 



For the semiclassical dynamic results presented ahead, an implementation of the HK 



propagator 



541 ] is used. It is reviewed below as is the linearized semiclassical initial value 



representation (LSC-IVR) 68|, l69[, which served as our workhorse to produce purely classical 
results to compare with the semiclassical and exact quantum ones. 



A. The Herman-Kluk Propagator 



The HK propagator is frequently used in chemical and molecular physics when treating 

ly. Its origins go back to 



53| and as an initial value 



systems with a large number of degrees of freedom numerica 
the work of Heller on frozen Gaussian wave packet dynamics 
representation, no root searches for classical trajectories need to be performed, as is the case 
for the Van Vleck-Gutzwiller (VVG) propagator 2M, l29|]. This fact may be advantageous 
especially for systems with several degrees of freedom, where root searching may become 
cumbersome. Formally, however, the HK is closely related to the VVG propagator via a 



stationary phase approximation 



57|. 



For a (i-degree-of-freedom system, the HK propagator reads 



i^"^(x,t;x',0) 



{2TihY 



(x|^^(qi,pt)) ■ i?(qo,Po,t) e'^"' ((?^(qo,Po)|x). (12) 
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Its basic ingredients are coherent states given in the position representation 



2V^'^^-7(x-q)V2+*p-(x- 



(x|^?,(q, P)) = ( - ) ' e-'^--^^^r2+^r>■i^^^)/n^ (^3) 



TT 



with the d X d width parameter matrix 7, which we assume to be diagonal and equal 
for (x|(7^(qt, pt)) and {gj{c[o,po)\x.'). Classical dynamics enters the expression through the 
propagation of trajectories from (qo,Po) to {qt,Pt) using Hamilton's equations, and the 
calculation of the classical action S = C dt', with C being the Lagrangian. An important 
role is played by the prefactor 



^(qo, Po, ^) = ^ y det l^mii + - ^^ni2i + "~— 12^ ' ("^'^^ 
with the elements of the stability (or monodromy) matrix: 

/m,, ^ gf\ .^g. 

For Gaussian wave packets, correlation functions can be computed directly with the HK 
propagator, according to 

c,,it) = I d'^xvl/*(x)vl/,(x,t) 

= 11 ^^^^0^ {^kito)\9,iqt,P^^^ (16) 
thus avoiding explicit calculation of the wave function. A 2(i-dimensional phase space integral 



remains which is usually calculated via a Monte-Carlo integration [70||. 

It is worthwhile to mention that the HK propagator is the lowest order term in a series 
expression for the full quantum mechanical propagator [tiI, [3] • Considering higher order 
terms in this series, the description of (deep) barrier tunneling is within the reach of semi- 



classical IVR methods 



73I . |7J] . The numerical effort to determine the correction terms is 



quite formidable, however. 

B. The Linearized Semiclassical Initial Value Representation 

The LSC-IVR serves to produce classical results and is summarized very briefly. It was 
developed by Miller and co-workers as an approximation to semiclassical IVRs, avoiding 
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oscillating integrands 68|, |69|]. Starting with an expression for the transition probability, 



i.e. the squared absolute value of the correlation function, utilizing a van-Vleck-type IVR 
and linearizing the phase differences of the trajectories one finally ends up with 



75|, 



where in this case At^) is the Wigner-Weyl transform of the associated initial wave packet 



^M(qo,po) 



-(«po-q)/?i 



q 



(qo + 77l*.-)(^j-|qo-77)- 



(18) 



Similarly, i?(^)(qf,pt) is defined using the final wave packet. Inserting explicit expressions 
for the wave packets \1/ gives 



Pik^-t) = (^) j I d"god%exp 



;^(Po-Pi; 

- 7(qo - q,)^ 



2vPi-Pfe)^-7(qt-qfc)^ 



(19) 



which is much easier to apply than the HK method. However, all phase information of the 
trajectories is lost and therefore Eq. (fT9|) represents a purely classical result. 



IV. THE MODEL SYSTEM 

As a model system for calculations, consider the two-dimensional quartic oscillator whose 



Hamiltonian is given by 18| 



i7(p, q) = ^i±^ + « (I + + '^Mhi) , (20) 

where A is the coupling coefficient and a = 0.2 a scaling factor. The factor b allows one 
to tune the symmetry of the system. We use b = it /A to have the convenient rectangular 
symmetry and secondly to be sufficiently far away from a "stripe" symmetry, which would 
involve a quasi one-dimensional system. This system is quite suitable because it reveals all 
kinds of dynamical behavior depending on the value of A. It can be tuned continuously from 
the fully integrable case via the mixed case towards the fully chaotic case in the sense, that 
nearly all trajectories explore the entire phase space for arbitrarily long times. 

For the visualization of a 4-dimensional phase space, a Poincare section (or surface of 
section) is a standard and convenient tool. For each of the three dynamical cases shown in 
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FIG. 1: (g'i,pi)-Poincare sections illustrating the chaotic domain for the coupling coefficients A 



0.0, A = —0.35 and A = —0.6 from top to bottom respectively. Reprinted from 



181 ] with permission. 
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FIG. 2: {q2,P2) and ((/i,pi)-Poincare sections for A = —0.35 showing the different subregions in 



the chaotic sea and t 
A; modified from 



le regular islands; the large regular island inside region 5' is denoted as island 



la] and reprinted with permission. 



Fig. [H the (gi, pi)-Poincare sections {p2 > 0) are for a single trajectory. The first case shows 
regular dynamics and corresponds to the situation of two uncoupled quartic oscillators. The 
regular motion (for other individual invariant tori) would appear as concentric curves around 
the origin. As A is decreased, the dynamics become mixed with near-integrable islands 
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(whose dynamics are governed by the KAM-theorem) surrounded by a chaotic sea. Arriving 
at A = —0.6, the whole phase space can be regarded as being fully chaotic (99.5% of the phase 
space is chaotic). A further decrease to A = —1.0 does not change the dynamics dramatically 
with the volume of the chaotic region being virtually identical with the phase space volume. 
This is also the borderline case which separates the bounded from the unbounded system 
and it is an interesting example of a system with a non-compact phase space, but one which 



nevertheless has only a discrete spectrum 



76|. 



The case of greatest interest in this work is the mixed case with A = —0.35. For this 
system, the various chaotic subregions separated by transport barriers formed by segments 
of certain stable and unstable manifolds (as well as the regular islands) were identified earlier 



m 



and the relevant figure is reproduced in Fig. [21 Classically, there are eleven subregions, 
however four of them are related by a reflection symmetry. Building the symmetry into the 
quantum mechanics effectively leaves only seven distinct subregions (related regions are 
denoted by a prime, as in 5 and 5'). The various phase space volumes and connecting fluxes 
are all known. As a final note, ahead time is measured in periods of the central closed orbit 
of the island A at the mean energy of the wave packets being used. Roughly then, t = 10 
can be thought of as the time a trajectory crosses back and forth across the width of the 
system a full ten times. 



V. CLASSICAL, SEMICLASSICAL, AND QUANTUM RESULTS 

Although the primary goal is to investigate whether HK implementation of semiclassical 
theory is sufficient to deduce localization information, it is valuable to begin by considering 
a case in which a classical transport barrier has no crossing flux. Amusingly, this may lead 
to the opposite situation of transport suppression. In this case, quantum propagation may 
reach classically forbidden subregions by way of tunneling. It turns out that the zero-flux 
regime can provide at least a couple of critical pieces of information. First of all, for quan- 
tum transport suppression to be understood properly, it is necessary to know how well a 
wave packet can be contained within a subregion. A wave packet does have exponentially 
suppressed tails that extend beyond whatever region is being considered. Ideally, their total 
weight external to their primary subregion should be on a far, far smaller scale than the 
extent of transport suppression one is trying to observe. One would thus know in which 
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regimes those tails are safely ignored. Cross-correlation functions evaluated to intermediate 
propagation times in a zero-flux regime reflect the measure of the tails, that is at least up to 
times where tunneling begins to take over. Secondly, the HK semiclassical propagator can 
seriously and erroneously magnify highly unstable trajectories leading to inaccurate propa- 
gation in strongly chaotic dynamical situations. This problem tends to grow exponentially 
with propagation time in chaotic systems. The zero-fiux regime allows one to see rather 
precisely at which times a particular HK propagation implementation loses accuracy. There 
the trajectories associated with the tails abruptly become magnified orders of magnitude 
beyond their true contributions, and the onset is readily seen. Increasing the number of 
trajectories in the implementation should then be seen to delay the onset of errors. 



A. The zero-flux regime 

The only zero-flux barriers in the quartic oscillators' dynamics effectively are the bound- 
aries of the regular islands, but they can be used for our purposes. The idea is to propagate a 
wave packet initially as well localized as possible within one such island and check its overlap 
with a wave packet localized as well as possible in its neighboring region. We choose the 
initial Gaussian wave packet inside island A and measure its overlap with a final Gaussian 
wave packet within region 5'. 

It is known that for this quartic oscillator system, chaos-assisted tunneling Q], which is 
a specific mechanism of dynamical tunneling 5|, is present for all the islands. This leads 
to an increasing probability for the propagated quantum state to have an overlap in the 
chaotic sea, but on a very long time scale, much longer than the local mixing time. Since 
the HK propagation is not implemented with tunneling trajectories, it should not show 
this increase at all, should be less than the quantum probability, and should be close to 
the classical result. The classical, HK semiclassical, and quantum magnitudes of the cross 
correlation function are plotted in Fig. [31 It is worth pointing out several features. Of 
primary importance is that the scale of the cross- correlation function is 0(10"''), which is 
much more than small enough for the purposes of studying suppression of transport ahead. 
It is also possible to see the increasing overlap with time due to the tunneling. Next it is 
worth commenting on the HK implementation. Considering the geometry of the semiclassical 
calculations, it is sufficient to choose a phase space annulus near the boundary of the regular 
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FIG. 3: Classical (solid horizontal line), semiclassical (long dashes) and quantum mechanical (dot- 
ted line) cross-correlation functions for the regular island with a probe gaussian in the chaotic 
region for different numbers of trajectories (tr.) in the HK calculations: 2 • 10^ tr. in (a), 64 • 10^ 
tr. in (b). The classical calculations were performed with 10^ trajectories. The quantum result 
was obtained with the split-operator fast Fourier transform technique [ttI . 



island mainly in the chaotic region for the initial conditions of trajectories in Eq. f[T^ . The 
contribution from the inner part of the wave packet is vanishing as those trajectories remain 
trapped well within the island. The upper limit of the ring is taken in such a way that the 
contribution from all integration points beyond this limit is negligible. Nevertheless, even 
with this simplification, 2 x 10^ trajectories in the HK calculation is not yet converged to 
below the tunneling contributions at the intermediate time scale shown. This behavior is 
cured by going to 64 x 10^ trajectories. This is a very stringent test in a strongly chaotic 
system and again well beyond what is necessary for the transport suppression study. To be 
a little more specific about the problems that chaotic systems create, when a part of the 
integration region is chaotic, two features work together to drastically hinder the convergence 
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of the semiclassical integral. Both originate in the sensitive dependence of the trajectories 
on their initial conditions. First of all, a set of initial points limited to a small part of 
the phase space still spread over the entire phase space. This implies a highly oscillatory 
integrand. Consequently one needs many more initial points than in the integrable case to 
obtain a sufficient coverage of the final phase space. Secondly, the prefactor in Eq. (1121) is 
a measure for the stability of the trajectories and it becomes exponentially large in time. 
It is the contributions of a locally dense set of trajectories with various phases that lead to 
cancellation of most of this magnitude. Again, it requires a number of sampling points being 
of the same order of magnitude to preserve the norm. Therefore, to improve the convergence 
of the integral over the annulus we employed the time-integration method presented in 78|. 



Also a cutoff value for the prefactor was used, as proposed in 
10% of the trajectories were discarded after a certain time. 



78 



79| . In this way about 



B. Partial transport barriers 

Moving on to the transport with non- vanishing fluxes, the time-averaged transition prob- 
ability is the physical quantity of interest. Several different examples are calculated that 
illustrate crossing, one, two or more barriers. In each case, the initial wave packet is the 
same, in the region 5' as indicated in Fig. [2], and the mean energy is in the neighborhood of 
the 1000th excited energy eigenstate. It turns out to be extremely convenient for comparison 
purposes to arrange that the smoothed local spectra of all the final wave packets considered 
have the same shape and width. In this way, the integral over the local spectra in Eq. ([6]) is 
identical in all cases. The criterion needed to ensure this is that the short time decay of the 
corresponding autocorrelation functions, in particular its absolute value, has to be equal. 
That simplifies the ratio of any two time-averaged transition probabilities to be 



(21) 



V,i A,, (A) h 

In order to arrange an appropriate initial decay behavior, it is sufficient to use the frozen 



Gaussian approximation 53|, assuming that on this time scale the width of the Gaussian 
changes only slightly. Expanding the potential up to first order, solving the canonical 
equations, and using the analytic result for Gaussian integrals with equal widths yields 



Pkk{t) ^ exp 



(22) 
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FIG. 4: Quantum mechanical transition probabilities in the zero barrier case (a) and for one barrier 
(b) at long propagation times. For the zero barrier case, both wave packets are in region 5' and for 
the one barrier case, the final wave packet is placed in region 6 (the regions are shown in Fig. [2]). 
The initial wave packet width is 75/ = 3.3. 



where Vq^ denotes the vector of the partial derivatives of the potential with respect to the 
positions at the initial point. As only the short time behavior is relevant, Eq. (122]) can be 
expanded to the quadratic time terms: 



(23) 



Assuming the mean energies of two wave packets are equal, then the smoothed local spectra 
cover the same energy range if Xj = Xk- Given the width of one wave packet, say the initial 
one 7~^, the width of the other, 7^^, is found by the relations 



7fc 



qo.fc 



(24) 
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if = and 

7^ = (v^. + ^^Jpo. + y(n. + ^^K.y-4%'^X.-p^,.) (25) 

otherwise. 

The transition probabihty for a case with no barrier crossing can be arranged by placing 
the final wave packet into the same region as the initial wave packet such that their initial 
overlap is negligible. To arrange a one barrier crossing example, the final wave packet is 
placed into the appropriate neighboring region. This leads to results of the kind shown 
in Fig. m There is a visible decrease in the locally averaged transition probability for the 
one-barrier case relative to the zero-barrier case. 

For a clearer visualization of the information about transport suppression, the transi- 
tion probability is smoothed with a Gaussian window function [as opposed to the "square" 
window used in Eq. ([3])] given by 

^..(to; = ^ / dte-(*-*o)^/(^^^)p,.(t), (26) 

— oo 

where C is a smoothing parameter. The value of C, is ideally chosen to be as small (narrow) as 
possible while still being capable of smearing out the rapid transition probability oscillations. 
Although the smoothing is not infinitely wide, it nevertheless should approach the ratio given 
in Eq. (|2T|). 

In Fig. [5], the quantum and classical time-smoothed transition probabilities are compared 
for the zero-barrier and one barrier cases. The shorter time dynamics in the upper panel 
illustrates the transient behavior and the bottom panel illustrates the long time average 
properties. In the transient behavior, the effect of the barrier is to delay the onset of 
the transition probability for both quantum and classical behaviors. In the theory given 
the one-barrier quantum and classical behaviors are both predicted to have a linear 



m 



increase initially, albeit with superposed oscillations. It is this time regime in which the 
HK propagation needs to be accurate if it is to be able to contain the information about 
long-time localization. At long times, the classical results cannot show any localization 
distinction and the ratio of the value of those two curves must approach unity by the time 
that mixing becomes complete. One sees clearly that this is not what happens at long times 
for the quantum curves and that the localizing effects due to the barrier prevents their ratio 
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FIG. 5: Locally time-smoothed quantum and classical transition probabilities for zero (j = 5' ,k = 
5') and one barrier (j = 5' ,k = 6). The quantum zero-barrier case is represented by the solid line, 
the classical zero-barrier case by the dashed line, the quantum one barrier case by dash-dots, and 
finally the classical one-barrier case by the dotted line. The smoothing parameter ^ = 0.35 in (a) 
and 0.8 in (b) and the results shown are obtained with 10^ trajectories. Note that the two classical 
cases are very close to each other at longer times whereas the two quantum cases are not. Also 
note the logarithmic scale in (a). 



from approaching unity. Crossing more barriers only increases these effects. In Fig. [HI a 
wave packet is placed in region 2 [see Fig. [2], which means that three barriers are effectively 
being crossed, and the transition probabilities are compared with the zero-barrier case shown 
in the previous figure. At short times, the onset delay is much clearer, but nevertheless, 
remarkably similar in the quantum and classical cases. Here, according to the approach 
of 60|], the initial increase is expected to be cubic in time. This follows by the tri-diagonal 
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FIG. 6: Locally time-smoothed quantum and classical transition probabilities for one (j = 5',k = 6) 
and three barriers {j = 5',k = 2). The quantum one-barrier case is represented by the solid line, 
the classical one-barrier case by the dashed line, the quantum three-barrier case by dash-dots, and 
finally the classical three-barrier case by the dotted line. The smoothing parameters and numbers 
of trajectories are identical to Fig. [SI Here the two classical cases take much longer to approach 
each other, but eventually do. On the other hand, the ratio of the values of the two quantum cases 
are far from unity indicating the extra localization due to three barriers rather than one. 

nature of the flux matrix, A, for this case and the expansion 

P(t) = exp(-At)P(0) = (^1 - At + ^AH^ - ^AH^^ P(0) + 0{t^) (27) 

Until the A matrix is cubed, the element of exp(— At) connecting one region to another 
across three barriers is null. In the long time dynamics, again the two classical transition 
probabilities approach each other, but it does take more time, and that is just as expected. 
The ratio of the local mean values of the quantum curves is again far from unity, which 
implies further localization also as expected. 
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FIG. 7: Comparison of the semiclassically and quantum mechanically calculated transition proba- 
bility. Shown is the one-barrier case in the previous two figures: quantum result (solid line) versus 
semiclassical result with 10^ trajectories (dashed line) and semiclassical result with 10^ trajectories 
(dotted line). The convergence for longer times with greater numbers of trajectories is clear. 

Next, consider the time- dependent semiclassical description of the zero and one barrier 
dynamics. First of all, it is important to understand the convergence behavior of the HK 
propagator given the presence of strongly chaotic dynamics. As an example, the semiclassical 
transition probability for the one-barrier case is compared with the quantum transition 
probability for different numbers of sampling points in Fig. [3 The numerical integration 
was performed with the importance sampling Monte-Carlo method 80|] . Note that the time 
integration method was not applied here, but a cutoff for the prefactor was used, too. It leads 
to discarding around 15% of the trajectories at the final integration time. Up to time t = 4, 
there is no significant difference between the semiclassical curves for a different number of 
sampling points. After that time, it is necessary to use the larger number of trajectories, 
10^, to achieve reasonably close convergence to the quantum result, although a complete 
convergence is not yet attained. 

The convergence with the greater number of trajectories turns out to be sufficient for the 
demonstration of localization due to transport barriers, but having sufficient trajectories is 
absolutely critical. In Fig. [6l it is seen that the transient rise to the long time average occurs 
mostly in the time range 1 < t < 4 for the one-barrier case and 3 < t < 7 for the three- 
barrier case. It is this time range from which the long time localization can be deduced. 
The greater number of trajectories does give convergence there as seen also in Fig. [8] for the 
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FIG. 8: Comparison of the semiclassical and quantum mechanically calculated transition proba- 
bility for the three barrier case: quantum result (solid line) versus the semiclassical result with 10^ 
trajectories (dashed line) and the semiclassical result with lO"^ trajectories (dotted line). 

three barrier case. However, a simple renormalization of the wave function for low numbers 
of trajectories as done in 40|| would not incorporate this behavior and would be insufficient 
for understanding localization. 

We now arrive at the crux of the matter, i.e. whether a realizable implementation of 
the HK propagator can contain the information necessary to deduce localization properties 
and whether that information can be extracted. To understand the results for the one- 
barrier case quantitatively, consider first the ratios of the long-time averages of the transition 
probabilities in the quantum and classical cases. The exact classical ratio always approaches 
unity; in our numerical results (see Fig. [5]) a close value of Pse/'Pss = 0.89 is found. In 
the next paragraph, the many trade-offs made in order to set up the test given here is 
addressed and how that can lead to the classical result not being exactly unity. Given this 
classical result, roughly speaking, a 10% variation due to all the competing conditions is 
about the best accuracy that can be expected. Back to the quantum result, which reflects 
the localization of interest: the further from unity, the greater the localization. For the 
one-barrier to zero-barrier ratio, the quantum calculations give for the ratio 0.34, which 
is clearly well away from unity to the level of variation that can be controlled, and it 
reflects the effects of the transport barrier. The next step is to flnd out how this compares 
to the semiclassical predictions. There are two ways to do this, and they only rely on the 
intermediate time dynamics to work. The flrst possibility is to use the perturbative analytical 
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results incorporating the measured classical fluxes, relative volumes of the subregions, and 
the mean excitation energy (in order to obtain the smoothed level density), which can all 



be obtained from [18|, in Eq. (PT]) together with Eqs. (ITIS]) . Note that if the system is 
not in the perturbative regime, the Master equation approach would have to be used. The 
other way possible in the perturbative regime includes the usage of the results from the HK 
propagation. We have shown that the HK method is capable of reproducing the behavior 
in the transient regime correctly, although great care must be exercised in order to ensure 
that enough trajectories are included in the implementation. To obtain information about 
the long time behavior, one can deduce numerically an approximate slope from this regime 
(see again Fig. 7), which corresponds to the linear coefficient of the transition probability 
in Eq. This again can be used to extract A56 respectively A56. Together with A54 from 
the perturbative analytical results one obtains A55 and furthermore the ratio. These two 
methods give 0.38 and 0.30, respectively, and bracket the actual quantum result. They are 
within roughly 10% of the quantum result. To the level that could have been expected under 
the conditions of the calculation, it indeed was possible to extract the information from the 
results using the HK propagator. 

Finally, it is worth highlighting explicitly and all together in one location, the nature 
of a number of the desires and approximations required to obtain the results in the above 
paragraph. It is desirable to remain on an energy surface that is excited sufficiently that 
semiclassical approximations are valid. In particular, in addition to the usual semiclassical 
considerations, localizing a wave packet within a particular region in this problem improves 
as the ratio E^^'^/h increases. The smaller a region or the more filamentary, curved, or 
possessing of small regular islands, the more difficult it is to locate a wave packet wholly 
within that region and have the tails avoid regular islands. This argues for increasing the 
energy or conversely shrinking h. Working in the opposite direction is the desire that the 
example not be so excited as to be way beyond the regime that might be of interest in physical 
chemistry. That desire also has the added benefit of keeping the transition parameters within 
the perturbative regime for the particular quartic oscillators studied (A = 0.35). Much higher 
in energy than the example given here and that would no longer have been possible. Making 
the entire enterprise more complicated is the use of wave packet widths such that their 
strength functions had the same energy spread in order to simplify the interpretation of the 
final results. The end result is that the wave packet tails did cross subregional boundaries. 
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including regular islands. This explains the classical ratio above not approaching unity, 
but rather 0.89. Similarly, competing effects are going on in the quantum case where the 
agreement is quite good, but at least part of the quality of the agreement is likely due 
to canceling opposing effects of wave packet tails crossing boundaries between two chaotic 
subregions (inhibiting a finding of localization) and the effects of tails crossing into a regular 
region (enhancing a finding of localizing effects). 



VI. CONCLUSIONS AND OUTLOOK 



The principal concern in this paper is whether an initial value representation of the 
semiclassical propagator heavily relied upon in molecular and atomic contexts, the Herman- 
Kluk propagator, can be implemented well enough to capture the information necessary 
to understand subtle long-time quantum transport effects, such as wave packet localization. 
The "freezing out" of wave packet exploration or transport of energy is potentially extremely 
important. In this quest, for many reasons the possibility of success can only arise if prior 
analysis allows one to deduce the long-time effects from much shorter time dynamics. Not 
only does the Herman-Kluk propagator fundamentally lose accuracy with increasing time 
scales, an implementation goes more and more beyond practical computational limits, and 
in addition, is made more difficult as a system's dynamics becomes more strongly chaotic 



and as the number of degrees-of-freedom increases. As was shown previously 18 



the wave 



packet localization induced by classical transport barriers such as cantori or crossings of 
stable and unstable manifolds provides one arena in which the relatively short-time and 
long-time behaviors are intimately linked and understandable. 

A system of two coupled quartic oscillators was investigated as a simple, convenient, 
previously-explored example of coupled anharmonic oscillators that could be expected to 
contain similar dynamical features found in models of molecular and atomics systems (al- 
though perhaps they are more strongly chaotic than typically found in these latter models 
and thus more challenging in that sense), and in particular, they display a variety of classi- 
cal transport barriers. Another helpful feature is the existence of islands of regular motion. 
They could be used effectively as zero-flux transport barriers, which in turn provided a 
rather exacting test of the HK propagator implementation. Although, dynamical (actually, 
chaos-assisted) tunneling is present, it is weak enough to be orders of magnitude smaller 
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than the effects of the classical transport barriers of interest. By leaving out tunneling tra- 
jectories in the HK implementation, convergence could be imposed to below the importance 
of the tunneling contributions seen in the full quantum dynamics. 

At least in the quartic oscillators, a fairly strongly chaotic two-degree-of-freedom sys- 
tem, it turned out to be possible to obtain accurate propagation throughout the short to 
intermediate times needed to deduce the long-time wave packet localization. The amount of 
localization deduced from the HK calculation matched the actual quantum result to within 
10%. Perhaps, greater implementation efficiency and less unstable dynamics would allow 
for the treatment of systems with a few more degrees-of-freedom. It would be worthwhile 
to understand where the upper bound in treating more degrees-of-freedom is. 

As a final point, there has been a Wigner function based study of recurrences using 
the PuUen-Edmonds Hamiltonian, and it has revealed the usefuUness of classical trajectory 
approaches in another mixed phase space system [sil. Therefore, there is every reason to 
expect that it should be possible to study localization semiclassically in non-energy-scaling, 
a-bit-more-realistic systems. 
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